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> ■ 1 Introduction 
(N 

' The celebrated Di Perna-Lions theory, which was initiated in [7], claims that the well-posedness 

'sj" . of the transport equation implies the existence and uniqueness of the quasi-invariant flow gen- 

! erated by a Sobolev vector field, provided its divergence is bounded. See [5] for the case of 

' Gaussian measure as the reference measure. This theory has subsequently been extended to the 

case of BV vector fields by Ambrosio [1, 2], to the infinite dimensional Wiener space in [3, 9], 
and to the case of SDE in [16, 10], inspired by the work of Crippa and de Lellis [6]. In two recent 
papers [12, 13], Le Bris and Lions used again the ideas of the Di Perna-Lions theory to study 
rS ■ the Fokker-Planck type equations, showing the existence and uniqueness of solutions when the 

. coefficients have a certain Sobolev regularity (see also [11, 15] for some related results). In [14], 

the author studied the well-posedness of Fokker-Planck type equations on the Wiener space, 
mainly under the exponential integrability of the gradients and divergences of the coefficients. 
For the study of Fokker-Planck equations in the Hilbert space, see [4] and the references therein. 
The present work is motivated by [13] . We consider the following Fokker-Planck equations 

dtp + di{pbi) - dfj{aijp) = 0, (1.1) 

where the coefficients h = [hi,-- - ,bn) and a = {aij)i<ij<n have only weak regularity on the 
spatial variables, e.g. Sobolev or BV regularity. For the motivation of studying this kind of 
equations with irregular coefficients, see [13, Section 4]. The above equation (1.1) is closely 
related to the stochastic differential equation 

dXt = a{Xt)dWt + b{Xt)dt, Xo = x, 
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in which a = {(yik)i<i<n,i<k<in is a matrix- valued function such that a = \(Ta* (here a* is the 
transpose of o"), and Wt is an m-dimensional standard Brownian motion. In the following we 
always assume that the coefficient a has such a form, as in [13]. The adjoint equation of (1.1), 
called the backward Kolmogorov equation, reads as 

dtp - hdip - aijdfjp = 0. (1.2) 

We will mainly focus on the Fokker-Planck equation of divergence form (a notion introduced 
in (5.8) of [13]): 

dtp + di{pbi) - -di{aik(Tjkdjp) = 0. (1.3) 

As pointed out in [13, Subsection 5.1], this equation is relatively easier to deal with, compared to 
the Fokker-Planck equation (1.1) and the backward Kolmogorov equation (1.2), since its second 
order term di{aik(Tjkdjp) has a self-adjoint form. To state the main result of this paper, we 
introduce some notations. For a locally integrable function / : R, we say that / has 

BV Iqc regularity, denoted by / € BV;oc(R"), if for every i — 1, • • • ,n, the distributional partial 
derivative dif of / is a locally finite Radon measure on R". /~ := —(/AO) denotes the negative 
part of the function /. For a BV^oc vector field b : R" — ?> R", we write D ■ b = "^^^i dibi for its 
divergence, which is still a Radon measure. If Z) • 6 is absolutely continuous with respect to the 
Lebesgue measure dx, we denote by div(6) its Radon-Nikodym density. Our main result is the 
following 

Theorem 1.1. Assume that b and a satisfy: 

6€ (Li([0,T],BV,,,(R-)))", (Li([0,r],Li+L-(R")))", ^^^^^ 



div(6) € Li([0,r],LL(M'")), [div(6)]- G Li([0,T],L-(R-)); 



and 



a G (L2([0,r],Vr/;^'(R")))"'''", G (l2([0, T], + L~(R")))"^"'. (1.5) 

Then for each initial condition in n L°° (resp. fl L°°), the equation (1.3) has a unique 
solution in the space 

{pGL°°([0,r],L2nL°°) {resp. L°°([0, T], n L°°)), a*Vp e L\[0,T], L^)}. 

This theorem gives a positive answer to the open question raised in [13, Remark 12] on page 
1299. It also generalizes [11, Theorem 1.4], since the diffusion coefficient a in [11] is assumed 
to be independent on the spatial variables. The proof of Theorem 1.1 will be provided in the 
next section. As in [13, Section 7], we can adapt the above result to the other Fokker-Planck 
type equations (1.1) and (1.2), by imposing suitable conditions on 6^^ = 6 — ^div((70"*). Finally 
we give a brief discussion on the reason why we cannot further generalize it to the case where a 
has only BV regularity, see Remark 2.6. 

2 Proof of Theorem 1.1 

In this section, we present the proof of Theorem 1.1. We consider the Fokker-Planck equation 
of divergence form 

dtp + di{pbi) - -di{aik(Tjkdjp) = 0. (2.1) 
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First we give the mathematical meaning of the above equation. For a given initial condition 
po G L2 nL~ (resp. n L°°), a function p G L°°([0, T], n L~) (resp. L°°([0, T], n L°°)) 
satisfying a*Vp G L'^{{0,T], L'^) is called a weak solution to (2.1) if for all G C~([0,r) x M"), 
it holds 



npdt(pdxdt+ [ po(p{0,-)dx = - [ [ p{b,Vip)dxdt + ^ [ [ {a*Vp,a*V(p) dxdt. 

The existence of solutions to the above equation is the easier part, see for instance the 
beginning of [13, Subsection 5.4] for the case po £ (1 L°° , and [13, Subsection 6.1.2] for the 
case pq (z r\ L°° . Hence in the following we focus on the uniqueness part of Theorem 1.1 and 
follow the ideas in [13, Subsection 5.3]. The main difference is that a single convolution kernel 
is not enough to achieve our result, instead, we need a family of kernels as in the proof of [1, 
Theorem 3.5], see also [2, Theorem 5.1]. 

Let 

/C = jp G C7^°°(R",]R+) : supp(p) C B{1), I pdx = ll 

be the family of candidate convolution kernels, where B{1) is the open unit ball centered at the 
origin 0. Take p £ JC and define = e~'^ p{e~^ •) for e > 0. We regularize the equation (2.1) in 
the spatial variables: 

dtiPe *p) + Ps* di{pbi) - ^ps * di{aik(Tjkdjp) = 0. (2.2) 
We denote hy Pe = Pe * V and introduce the notation 

[Pe.c]{f) = p^ * {cf) - c{pe * f) 

for a differential operator c. Note that c can also be a real valued function. Using this notation, 
we have 

Pe * di{phi) = Pe* (div(6)p) + pe* [hdip) 

= [/Oe,div(6)](p) + [pe,bidi]{p) + di{hiPe) 

= Ql,e + Q2,e + d^{b,Pe), (2.3) 

where we have defined 

Qi,s = [ps,dw{b)]{p) and Q2,s = [pe,bid^]{p). (2.4) 

In fact, it is the term (52, e that causes the trouble in the BV situation and marks the difference 
between the present work and [13]. 
Next we have 

Pe * di{aik(7jkdjp) = Pe* {{dicrik)ajkdjp + crikdi{ajkdjp)) 

= [pe, diaik]{crjkdjp) + [diGik) Pe * {crjkdjp) 
+ [pe,crikdi]{ajkdjp) + aikdi{pe * (ajkdjp)). 

Define 

Se = [pe,diaik]{crjkdjp) and = [pe,crikdi]{ajkdjp). (2.5) 

We remark that the term Ss prevents us from extending the results to the case where the 
diffusion coefficient a has only BV regularity, see Remark 2.6 for more details. Now we have 

Pe * di{aik(7jkdjp) = Se+Ts + di{aikPe * {(Tjkdjp)). 
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If we denote by 

Rk,e = [Pe,crjkdj]{p), k = l,---,m; (2.6) 

then 

Pe * {(^jkdjp) = Rk,e + (^jkdjPe- 

Summing up the above discussions, we arrive at the equahty (5.17) in [13]: 

Pe * di{aik(7jkdjp) = Se + Te + di{aikRk,e) + di{aik(7jkdjPe). (2.7) 

Combining (2.2), (2.3) and (2.7), we get an equation of p^ similar to (2.1), but with some error 
terms on the right hand side: 

dtPs + di{pebi) - ^diiaikCTjkdjPs) = -Ql,e - Q2,e + ^{diiaikRk,e) + Se + Te). (2.8) 

Now we need the classical commutator estimate in the Di Perna-Lions theory (see [7, Lemma 
II. 1] or [13, Lemma 1]). We include it here for the sake of the readers' convenience. 

Lemma 2.1 (Commutator estimate I). For r, a, ri, ai > 1, set ^ = + ^ and ^ = ^ + -^- 

Let f € L'-i([0,T],LL,(R")), g e ([0, T], Lf,,(M")) and c G (L^^ ([0, T], H^/„'^"(M")))".' Then 
as e ^ 0, 

[Pe,cA]{f)^0 ^nL^l([0,^],Lf,^(]R")), (2.9) 

and 

[p„g]{f)^0 znL^^{[0,T],Ll^{W^)). (2.10) 

Remark 2.2. We observe that under the assumptions of the above lemma, for any fixed e > 0, 
the commutator [p^, Cidi]{f) belongs to the space L^i ([0, T],WI^^{W)) , i.e., it has the first order 
Sobolev regularity with respect to the spatial variables. This can be seen from its expression 
(see (3.8) in [1] or [8, Lemma 2.5]): 

[p„ Cid,]{f){x) = -[(/div(c)) * pe]{x) + / f{y){c{y) - c{x), {Vpe){x - y)) dy 

= -[(/div(c)) * pe]{x) + [{fc) * {Vpe)]{x) -{Cj* {Vpe)){x). 

Indeed, the three terms (/div(c)) * ps, {fc) * (Vpe) and / * (Vpe) are smooth, therefore the 
assertion follows from the fact that c G (L°i ([0, T], VF/;^''(M")))". □ 

We concentrate in the following on the L^-theory, that is, we prove the uniqueness of solutions 
to (2.1) in the space 

X2 = {pe ([0, T] , l2 n L°°) ; a* Vp G ( [0, T] , L^) } . 

For the L^-theory, the same argument as in [13, Subsection 6.1.2] works. 

Now we consider the error terms in (2.8). Using Lemma 2.1, we see that for the time- 
dependent vector field b G (^^([0, T], BV;oc(]R")))" satisfying div(6) G ([0, T], ^^^^(M")) , the 
first error term 

Qi,, = [pe,divmp) '-^ in L\[0,T],lUr-)), (2.11) 

since p G L°° ([0, T], L°°(M")) . 

The estimate of the commutator Q2,e cannot be obtained from Lemma 2.1, instead, we will 
rely on the work of Ambrosio (see [1, Theorem 3.2] or [2, Section 5]). We first introduce some 
notations. Let Dbf = {djbi{t))i<ij<n be the "Jacobi" matrix of distributional derivatives of bt, 
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whose entries are locally finite Radon measures since bf E BV/oc(I^")- Denote by \Dbf\ its total 
variation. Let 

Dbt = D%t + D'bt 

be the Lebesgue decomposition of Dbt into absolutely continuous and singular part with respect 
to the Lebesgue measure dx. Define the matrix valued function M : [0,T] x (8) M" as 

the Radon-Nikodym derivative of Dh with respect to \Db\. Here the measure \Db\ on [0, T] x M" 
is defined as 



ip{t,x)(i\Db\{t,x) 



[ [ (p{t,x)d\Dbt\{x)dt, for ah 93 G Cc([0,r] X M"). 
Jo Jm."- 



Db, {D^-bl and iD'^bl are defined similarly. We remark that Ambrosio defined in [1] the matrix 
Mt to be the Radon-Nikodym derivative of D^b with respect to |-D*6|; however, when restricted 
on the support of |i:>^6|, we have Mt{x) = Mj(x) for \D%\-a.e. {t,x) G [0,r] x M". 

Lemma 2.3 (Commutator estimate II). For any compact set K C (0, T) x M", we have 

m[ \Q2,e\dxdt<\\p\\^ [ A{Mt{x),p)d\D'b\{t,x) + \\p\U{n + I{p))\D%\{K) (2.12) 
£4^0 Jx Jk 



and 



lim / \Q2,e\dxdt<\\p\\^I{p)\D'b\{K), (2.13) 
where for a matrix M and p G C^(]R"), 

/ \z\-\Vp{z)\dz. 



A(Af,p) = / \{Mz,Vp{z))\dz, I{p) 

Now we turn to the error terms concerning the diffusion coefficient a. The arguments are 
similar to those in [13, Subsection 5.3]. Applying again Lemma 2.1, we have 

Rk,e = [Pe,<yjkdj]{p)'-^Q inL2([0,T],LL(M")), (2.14) 

due to the facts that a^k G ([0, T], W/;^^(M")) and p G ([0, T], L°°(M")) . Next, since both 
diUik and ajkdjP belong to L"^ {[Q,T],Ll^{W)) , 

Se = [pe, diaik]{ajkd,p) '-^ in L\[0, T],lI^{R^)) . (2.15) 
Finally, as aik G ([0, T], w/^f (M")) , i = 1, • ■ • , n, we have 

T, = [p„aik^^]{a,k^jp)'-^0 m L\[0,T], lI^{M.^)) . (2.16) 

We denote by 

Ue = -Ql,e + l{Se + T,y, 

then the estimates (2.11), (2.15) and (2.16) lead to 

f/^i^O inLi([0,T],LL(M")). (2.17) 
The equation (2.8) can be rewritten as follows: 

dtPe + di{pebi) - ^di{aikajkdjPe) = 11^- Q2,e + ^di{aikRk,6)- (2.18) 
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Notice that is smooth with respect to the spatial variable. By Remark 2.2, the commutator 
Rk^e € L^([0, T], W;^'^(M")) , which together with the regularity assumptions on the coefficients 
h and a tells us that dtPe G ([0' ^] ^ ^'') • Therefore Pe G W/^'J ( [0, T] x M") and we can apply 
the standard chain rule in Sobolev spaces. For j3 € C^(M), one has 



dtP{Pe) + di{P{Pe)h) - 7;di{crikCrjkdjf3{pe)) 



1 

2' 

1 



= (3'{Pe)dtPe + bil3'iPs)diPe + /3(pe)div(6) - -di{aikCrjk(3'{pe)djPe) 
= P'{Pe)dtPe + P'{pe)di{peh) + [P{pe) " div(6) - ^0 {pe)di{a,k(T jkdjPe) 

- ^{(^ikCrjkdjPe)(3" iPe)diPe 
= P'ipe) [dtPe + d,{pebi) - \di{a,k<y jkdjPe)] + (/3(Pe) " div(6) - ^/?"(p,)|a*Vpe|'. 

By (2.18), we obtain 

dtpipe) + dmPe)h) - \d,{w,kdj(3{p,)) - {f3{p,)-p,p'{p,))divib) + ^p"ip,)\a*Vpe\^ 

= P'{pe) [Ue - Q2,e + \d^{aikRk,e)\ ■ (2.19) 

In order to prove the uniqueness of solutions to (2.1), we shall use the technique of renor- 
malized solutions, a notion which was introduced by Di Perna and Lions in [7, Section II. 3] (see 
also [2, Definition 4.1] and [11, Definition 4.9]). 

Definition 2.4 (Renormalized solution) . iei cr : [0, T] x M" ^ R"* (g)R" and 6 : [0, T] x M" 

M*" he such that (i) a G ^^([0' ^] ^ ^'')' N div(6) G Lj^^ii^, T] x M"). We say that a solution 
p to (2.1) is a renormalized solution if for any /3 G C^(M), the following equation holds in the 
distributional sense: 

dtm + ^^{mh) - \di{aika,kd,m) " (/5(p) " P/3'(p))div(6) + i/3"(p)|a*Vp|2 = 0. (2.20) 

Now we show that any weak solution in the space X2 of (2.1) is renormalizable, provided 
the conditions of Theorem 1.1 are satisfied. The main idea of the arguments are similar to the 
proof of [1, Theorem 3.5] (see in particular Step 3 therein). 

Theorem 2.5 (Renormalization property). Under the conditions of Theorem 1.1, any weak 
solution 

peX2 = {pe ([0, T], L2 n L~) ; a*Vp G ([0, T], L^) } 
of (2.1) is also a renormalized solution. 

Proof. We have to show that as e J, 0, all the terms on the left hand side of the equation (2.19) 
converge in the distributional sense to the corresponding ones in (2.20), while the limit of the 
right hand side is 0. We split the proof into two steps: in the first step we show the convergences 
of all the terms except the one involving Q2,e^ while in the second step we focus on the term 

0{Pe)Q2,e- 

Step 1. Since p is essentially bounded and /3 G C^(M), /3 and its derivatives are uniformly 
continuous on the interval Ip = [— ||p||oo, ||p||oo]- Notice also that ||pe||oo < IIpIIoo- As the proofs 
are similar, we only illustrate the convergences of 

diicrik(Tjkdj(3{pe)) di{aik(Tjkdj(3{p)) and (3'{pe)di{aikRk,e) ^ 
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as e J, in the distributional sense. 

First for any 93 € C^([0, T) x M"), we have by the integration by parts formula that 

/ / ifdi{aik(Tjkdj(3{ps))dxdt- / Lpdi{aikcrjkdjfi{p))Axdt 
Jo JR" Jo Jr" 

< r I \p'{pe){a*V^,a*Vpe) - P'{p){a*V^,a*Vp)\<ixdt. 

Jo JK" 

The above quantity is dominated by the sum of 

Jl,e ■■= ||/3'(Pe)||oo||V(^||oo / \a\ ■ \a*Vpe - a*Vp\ dxdt 

Jk 

and 

J2,e := ||V</.||oo / W{Pe) " {p)\ ' k| • \a*Vp\ dxdt, 



IK 

where the compact set K := supp(<y9) C [0,T] x M" is the support of (p. By Cauchy's inequality, 
we have 



1 , ,1 



Ji,e<c(^j \a\^dxdty (^j^ \a*Vpe-a*\/p\^dxdtY ^0 

as e 4, 0. In this paper C denotes the constants whose values have no importance and may 
change from line to line. Next since p G L°^{[0,T], L'^{W)) C L'^{[0,T], (W)), p^ tends to p 
in the latter space as e J, 0; consequently, (3'{pe) converges to P'{p) in measure. Moreover, 

\/3'{Pe) - /3'{p)\ • kl • \a*Vp\ < {\Wipe)\\oo + 11/3' (p) Hoc) kl • \a*Vp\ < C\a\ • \a*Vp\ 

which is integrable on K. Therefore the dominated convergence theorem tells us that 

lim J2,e = 0. 
£4,0 

To sum up, we conclude that when e 10, di{aik(Jjkdjf3{ps)) converges to di{aik(Jjkdjl3{p)) in the 
distributional sense. 

Now we consider the limit P' {pe)di{aikRk,e) ^ 0. For ip € C^[[0,T) x M"), again by 
integrating by parts, one has 

n<p (3' {pe ) di {(JikRk,e) dxdt 

= -11 f3'{p,){a*V^,Re)dxdt- [ [ ^p"{p,){a*Vpe,Re)dxdt 
Jo Jr" Jo Jr" 

where R^ = (i?i,e, • • • ,Rm,e)- Note that 

\J3,s\ < |l/3'(p,)||oo||V(^||oo / kl • \Rs\dxdt 

Jk 



<Ci \a\^dxdty (^j \R,\'^dxd?^ 



Therefore by (2.14), we get lim^^o "^a.e = 0- Since (3"{ps) and ||cr*Vpe||^2([o,T]xR") bounded, 
uniformly in e > 0, we can show that lim^j^o Ji,e = in the same way. 
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Summing up these discussions, we conclude that as e J, 0, all the terms in (2.19), except 
{Pfr)Q2^e-, converge in the distributional sense. We define the "defect" measure 

^x ■= ^t|3ip)+^^i|3ip)b^) - ^d,iaikCTjkd,l3ip)) - {l3{p) - p/3'(p))div(6) + ^f3"{p)\a*S/p\^ (2.21) 

(notice that j3'{pe) \Us + \di{aikRk,e)\ ^ as e tends to 0). 

Step 2. Now we deal with the term —/3'(p^)Q2^e- Let Qp be one of the weak limit points 
of {pe)Q2,e\ in the sense of measure (by Lemma 2.3, such an accumulating point exists since 
\P'iPe)Q2,e\ is bounded in L}J. The measure Qp depends on the convolution kernel p, but the 
"defect" measure p, defined in (2.21) is independent of p and satisfies \p\ < Qp. Thus we get 

\p\ < Qp for ah pelC. (2.22) 

We deduce from (2.13) that Qp is absolutely continuous with respect to which together 

with (2.12) gives us 

Qp<\\p\\^A{M.{-),p)\D'b\. 

Thus by (2.22) we obtain 

\p\ < \\p\\^A{M.{-),p)\D'b\, for all p G /C. 
Denote by g the Radon-Nikodym density of \p\ with respect to Then for each p G /C, 

g{t,x) < \\p\\ocA{Mt{x), p) for \D^b\-a.e. {t,x). 
Let D be a countable dense subset of IC with respect to the norm W^'^{B(1)). We have 

g{t,x) < IIpIIoo inf A{Mt{x),p) for \D'b\-a.e. {t,x). (2.23) 

For iD^^I-a.e. {t,x) fixed, we deduce from the definition of A{Mt{x),p) that the mapping 
IC 3 p ^ A{Mt{x),p) is continuous with respect to the W^'^{B{1)) norm. Therefore 

inf A{Mt{x),p) = inf A{Mt{x),p). 
pev p&K. 

Now by Alberti's rank one structure of Mt{x) (cf. [1, Theorem 2.3]), we conclude that the above 
infimum is for \D^b\-a.e. {t,x) G [0,T] x (see [1, Lemma 3.3]). Therefore 5 = \D^b\-a.e. 
by (2.23). As a result, the measure ^ = 0. By (2.21), we arrive at 

+ di{/3{p)h) - ldi{a,kCTjkdj/3{p)) " (/5(p) - P/3'(p))div(6) + ^l3"ip)\a*Vp\^ = 0. (2.24) 

The proof is complete. □ 

Finally we are ready to prove the uniqueness of solutions to (1.3). 

Proof of Theorem 1.1. Since the equation (1.3) is linear, it is enough to show that if p is a 
weak solution such that p{0) = 0, then p{t) = for all t >0. By Theorem 2.5, we know that for 
any (3 G C^(M), (2.24) holds in the distributional sense. If we choose /3(s) = for s G M, then 
it becomes 

dt{p^) + di{p\) - d^iaikajkpdjp) + /div(6) < 0. (2.25) 

We take (j) G C^(M'\ [0, 1]) satisfying (I)\b{i) = 1 and supp((/)) C B{2). Define a nonnegative 
smooth cut-off function ipR = 4>{-k) for R> 0; then \/(j)R = ^^<P{-r)- To simplify the notations. 
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we will write f for the integral of the function / on with respect to the Lebesgue measure. 
Multiplying the inequality (2.25) with (pji and integrating by parts on R", we obtain 



d 
dt 



/ /(6,V0r)-/' p{a*Vcl)n,(T*Vp)- [ pV/jdiv(6). (2.26) 

JR" JR" JR" JR" 



Now we estimate the three terms on the right hand side. First 

1 



{R<\x\<2R} 

<C^I|V<^||l 



R ^'R' 

2 1^1 



P 



{\x\>R} 



1 + Ixl 



Since G {L^{[0, T],L^ + L°°(M")))", there are two vector fields 6i, 62 such that b = bi + b; 

|i^eLi([o,r],Li(M")), 

Mt)\ 



and 4h e Li([0,r],Li(M")), G Li([0, T], L-(M")). Then we have 



jR" 



^ C'|[p|||oo([o^T],L°°) 



{|x|>fl} 1 + fI 



&2(t) 


1 


1 + \x\ 


LOO J{\ 



p 



2.27) 



Note that as a function of t G [0,T], the right hand side of (2.27) is dominated by 



C 



hit) 



1 + X 



+ C'IIpIIl°°([o,t],l2) 



hit) 



1 + X 



and the latter is an integrable function of t G [0, T]. Furthermore, for a.e. t G [0,T], the right 
hand side of (2.27) tends to as i? ^ 00. Therefore by the dominated convergence theorem, 



hm / / p\b,VcPR)=0 

Next, in a similar way. 



(2.28) 



{R<\x\<2R} 

<C^I|V0||l' 



|p|-k*Vp|--||V,/.||L^|a| 



k*vpj- 



<7 



'{\x\>R} l + \x\ 

By the assumptions on a, we can split it into two matrix- valued functions cri and CJ2, such that 



1 + |x, 
By Cauchy's inequality, 

/ p{a*Vp,a*V^R) 

JR" 



G L2([0,T],L2(]R")) and 



1 + \x\ 



G L2([0,r],L~(M")). 



<C\\p\\l 

< C\\a*Vp\\L2 
+ C 



{\x\>R} l + \x\ 



+ c 



F2 



{\x\>R} i + 



2 \ 1/2 



a2(t) 


/ 


1 + |x| 


L°° -'{1 



{\x\>R} (1 + 

bl • k*vp|. 



(2.29) 
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It is easy to see that the right hand side of (2.29) is dominated by 



C\\a*Vp\\L2 



1 + X 



L2 



+ C'l|p||L°°([0,T],L2)||cr*Vp||i2 



(12 (t) 



1 + \x\ 



which, by Cauchy's inequahty and the properties of o'i,o'2 and p, is integrable with respect to 
t € [0, T]. Moreover, for a.e. t € [0, T], the right hand side of (2.29) vanishes as R ^ oo. 
Consequently, the dominated convergence theorem leads to 



(2.30) 



lim / / p{a*Vp,a*V(t)R) = 0. 

Finally it is clear that 

p^cPRdw{b)< [ pVif [div(6)]- < ||[div(6)]-||Loo / (2.31) 
Now we integrate the equation (2.26) in time from to t (note that p(0) = 0) and obtain 



P 4>R< 



P^{b,V(l)R)- [ [ p{a*V(PR,a*Vp)- [ [ pViidiv(6). 

" JO JM" Jo JM" 



By (2.28) and (2.30), for any rj > there exists Rq > such that for any R > Rq, we have 



T 



p\b,V<PR) 



r [ p{a*V<PR,a*Vp) 

Jo JM" 



< 



Taking into consideration the inequality (2.31), we get 



/ /0i?<??+ r \\[diyib)]-\\L^ [ p\ 
Jr" Jo Jr"- 



First letting R t +oo and then r/ J, 0, we finally get 



p2 < 



T 



llV 



p 



Since [div(6)]- € ([0, T], L°°(]R")) , we have J^„p'^ = 0. This shows p = a.e. and the 
uniqueness is proved. □ 

Remark 2.6. Now we briefly discuss the reason why we are unable to deal with the diffusion 
coefficients a of BV regularity. Recall the definition of Ss in (2.5): 

Ss = [Pe,diaik\{crjkdjp) = Ps * {{crjkdjp)diaik) - {pe * {ajkdjp))diaik. 

For simplicity, we assume that the functions are time-independent and denote by g' = (Jjkdjp € 
L2(M"). Then 

Se = Pe* [gdiCTik) - {pe * g)diaik. 

If (Tjfc has only BV^oc regularity, then /x := SjCJjfc is a locally finite Radon measure on M". We 
have the decomposition ^ = V^fi + D^fi, where D'^fj, <C dx and D^fi ± dx. In the case D^fx ^ 0, 
since g G -L^(M") is not continuous, the product gD^fi is sensitive to the modification of g in 
Lebesgue negligible sets. Therefore is not a well-defined object on M". □ 
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